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Abstract. We study the renormalization group flow for the Higgs self cou- 
pling in the presence of gravitational correction terms. We show that the 
resulting equation is equivalent to a singular linear ODE, which has explicit 
solutions in terms of hypergeometric functions. We discuss the implications 
of this model with gravitational corrections on the Higgs mass estimates in 
particle physics models based on the spectral action functional. 



Ella estd en el horizonte. Me acerco dos pasos, ella se aleja dos pasos. Camino 
diez pasos y el horizonte se corre diez pasos mds alia. Por mucho que yo camine, 
nunca la alcanzare. ^Para que sirve la utopia? Para eso sirve: para caminar. 

(Eduardo Galeano) 

1. Introduction 

A realistic Higgs mass estimate of 126 GeV was obtained by Shaposhnikov and 
Wetterich in |35j . based on a renormalization group analysis, using the functional 
renormalization group equations (FRGE) method of [40 for gravity coupled to 
matter. In this setting, the renormalization group equations (RGE) for the matter 
sector acquire correction terms coming from the gravitational parameters, which 
are expressible as additional terms in the beta functions which depend on certain 
parameters: the anomalous dimensions a x and the scale dependence po of the New- 
ton constant. These gravitational correction terms to the RGE make the matter 
couplings asymptotically free, giving rise to a "Gaussian matter fixed point" (see 

E3: ESI). 

In this paper we carry out a more detailed mathematical analysis of the RGE 
for the Higgs self-coupling, in this asymptotic safety scenario with anomalous di- 
mensions, using the standard approximation that keeps only the dominant term in 
the Yukawa coupling matrices coming from the top quark Yukawa coupling. We 
show that the resulting equations have explicit solutions in closed form, which can 
be expressed in terms of hypergeometric functions. 

We discuss the implications of using this RGE flow on the particle physics models 
based on the spectral action functional. In particular, one can obtain in this way 
a realistic Higgs mass estimate without introducing any additional field content to 
the model (see the recent [H] for a different approach based on a coupling with 
a scalar field), but the fact that the RGE with anomalous dimensions lead to a 
"Gaussian matter fixed point" at high energies requires a reinterpretation of the 
geometric constraints at unification energy imposed by the geometry of the spectral 
action models, as in [IB] , 
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The paper is organized as follows. In £|2]we review the basic formalism of models 
of matter coupled to gravity based on noncommutative geometry and the spectral 
action functional, and in particular of the model described in [13] and in Chapter 
1 of [H]. In *J3] we review the use of renormalization group analysis in these mod- 
els, in the case where the renormalization group equations used are those of the 
Minimal Standard Model (MSM) or of the extension with right handed neutrinos 
with Majorana mass terms (z/MSM). We describe the usual approximations to the 
full system of equations, in view of our later use of the same approximations in the 
presence of gravitational corrections. We also show how it is natural to think, in 
these models, of the RG flow as a flow of the finite noncommutative geometry, or 
equivalently a flow on the moduli space of Dirac operators on the finite geometry. 
In Qwe show that it makes sense to apply Wetterich's functional renormalization 
group equations (FRGE) developed in [3D] to the spectral action. Using the high 
precision of the approximation of the spectral action by an action functional for the 
Higgs field non-minimally coupled to gravity, this leads to renormalization group 
equations of the type derived in [25) . |29j for matter coupled to gravity, within the 
asymptotic safety scenario, with a Gaussian matter fixed point. These renormaliza- 
tion group equations differ from the usual ones described in <|3]by the presence of 
gravitational correction terms with anomalous dimensions. In |5]we study the re- 
sulting RGE with anomalous dimensions, under the same approximations described 
in Sj3] We show that one can give explicit solutions for the resulting equations for 
the top Yukawa coupling and the Higgs self-coupling and that the latter is express- 
ible in terms of the Gauss hypergeometric function 2 Fi(a,b,c, z). We show that 
one can fit the Shapshnikov-Wetterich Higgs mass estimate within this setting. In 
^6] we discuss the effect of these modified RGE on the constraints imposed by the 
geometry of the model on the boundary conditions at unification. 



2. The spectral action and the NCG models of particle physics 

2.1. General geometric framework of NCG models. Noncommutative ge- 
ometry models for particle physics coupled to gravity are based on enriching the 
ordinary four-dimensional spacetime manifold X to a product X x F (more gen- 
erally, a nontrivial fibration) with "extra dimensions" F consisting of a noncom- 
mutative space. The noncommutative geometry F is described by a spectral triple, 
that is, data of the form (A,H,D), where A is an involutive algebra acting on a 
Hilbert space 'H, and a Dirac operator D on % satisfying a compatibility condi- 
tion with A, that commutators [D,a] are bounded operators, while D itself is a 
densely defined self-adjoint operator on % with compact resolvent. The spectral 
triple is additionally endowed with a Z/2Z-grading 7 on with [7,0] = and 
D7 = — 7-D, and a real structure J, which is an anti-linear isometry on % satisfy- 
ing J 2 = e, JD = e'DJ, and J7 = e"j J, where e, e', e" are signs ±1 that determine 
the KO-dimension of the noncommutative space. The real structure also satisfies 
compatibility conditions with the algebra representation and the Dirac operator, 
given by [a, b°] = for all a, b G A with b° — Jb* J -1 and the order one condition 
for the Dirac operator: [[D, a], b°] = for all a, b G A. 

In these models one generally assumes that the noncommutative space F is finite, 
which means that A and % are finite dimensional. In this case, the spectral triple 
data become just linear algebra data and the conditions listed above reduce to 
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equations D* = D, [a, b ] = 0, and [[D, a], b°] = 0, to be solved with the constraints 
[7, a] = 0, L>7 = -7D, J 2 =e, JD = e'DJ, and J7 = £"7 J. 

2.2. i/MSM coupled to gravity from NCG. In [T3] a particle physics model 
coupled to gravity is constructed with this method. It is shown (see also §1 of [16]) 
that it recovers an extension of the minimal Standard Model with right handed 
neutrinos and Majorana mass terms. The ansatz algebra for the finite geometry F, 
in the model of [T3] is taken to be Alr = C © H L © U R © M 3 (C). The main step of 
the construction of [T3| are the following. A representation of the algebra Alr is 
obtained by taking the sum M. p of all the inequivalent irreducible vAiij-bimodules 
(with an odd condition that refers to the action of a natural involution in the algebra, 
see [T3]). The number N of particle generations is not predicted by the model and 
is assigned by taking as Hilbert space of the finite spectral triple Tip — ® N Mf- 
This representation space provides the fermion fields content of the particle physics 
model. The left-right chirality symmetry of Alr. is broken spontaneously by the 
order one condition of the Dirac operator. This selects a maximal subalgebra on 
which the condition holds, while still allowing the Dirac operator to mix the matter 
and antimatter sectors. The subalgebra is of the form Ap — C © H © M 3 (C). The 
resulting noncommutative space F is metrically zero dimensional but with KO- 
dimension 6. The real structure involution Jp exchanges matter and antimatter, 
and the grading 7^ distinguishes left and right chirality of particles. There is a 
complete classification of all possible Dirac operators on this finite geometry and 
the parameters of the particle physics model that include Yukawa parameters Y 
(masses and mixing angles) and the Majorana mass terms M of the right handed 
neutrinos geometrically arise as coordinates on the moduli space of Dirac operators 
on the finite geometry F. The boson fields of the model arise as fluctuations of 
the Dirac operator D h-> Da, with the gauge bosons corresponding to fluctuation 
in the horizontal (manifold) directions and the Higgs field arising as fluctuations 
in the vertical (noncommutative) direction of the product space X x F. In recent 
supersymmetric versions of the NCG model [3], the supersymmetric partners of the 
Standard Model fermions also arise as vertical fluctuations of the Dirac operator, 
for a different choice of the finite geometry F. 

The field content of the model, including the neutrino sector, agrees with the 
z^MSM model discussed in [35] , [33J , [31] ■ This is an extension of the minimal Stan- 
dard Model (MSM) by right handed neutrinos with Majorana masses. In addition 
to the parameters of the MSM there are additional real parameters that correspond 
to Majorana neutrino masses, and additional Yukawa coupling parameters for the 
lepton sector given by Dirac neutrino masses, mixing angles, and CP-violating 
phases. The resulting model has a total number of 31 real parameters coming from 
the moduli space of Dirac operators on the finite geometry F (which account for 
all the parameters listed here) plus three coupling constants. A significant differ- 
ence with respect to the j/MSM model lies in the fact that this NCG model has a 
unification energy. 

2.3. The spectral action functional. The spectral action functional for the 
model of |13j is then given by 

(2.1) Tr(f(D A /A)) + ^{JlD A i), 
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where the first term gives rise, in the asymptotic expansion, to all the bosonic 
terms of the particle physics Lagrangian, coupled (non-minimally) to gravity, and 
the second term gives the fermionic terms, and their interactions with bosons. The 
variables £ are fermions in the representation %p (g> L 2 (X,S), with S the spinor 
bundle on X, viewed as anticommuting Grassmann variables. The asymptotic 
expansion for large A of the spectral action [10] is of the form 



(2.2) Tr(/(D/A))~ £ f k A k l\D\~ k + /(0)C D (0) + o(l), 

fcSDimSp+ 

where = J °° f(v)v k ~ 1 dv and /o = /(0) are the momenta of the test function /, 
and where the integration -ja\D\~ k is defined by the residues of the family of zeta 
functions Cd( s ) = Tr(|Z?|~ s ) at the positive points of the dimension spectrum of 
the spectral triple, that is, their set of poles. 



The asymptotic expansion (2.2 1, as computed in j!3j . gives terms of the form 
(2.3) 

Tr(/(D A /A)) ~ ^(48/ 4 A 4 -/ 2 A 2 c+^0) J ^gd'x 
96/ 2 A 2 -/ c 



^ J (~R*R* -3C W<T C^nV9d 



24?r 2 

10 

(-2a/ 2 A 2 + e/ 



A x 



12 



X 



+ ^2 / (53 2 ^ v + <? 2 F; v +\g\ B^ Bn ^gd 4 x, 



where R*R* is a topological non-dynamical term that integrates to a multiple of 
the Euler characteristic x(^), and C^ vpa is the Weyl curvature tensor of conformal 
gravity. The Higgs field is non-minimally coupled to gravity through the term -R|<y9| 2 . 
The first two terms provide cosmological and Einstcin-Hilbert gravitational terms 
and the remaining terms are the kinetic and self-interaction terms of the Higgs and 
the Yang-Mills terms of the gauge bosons. The coefficients of these various terms 
depend on the parameters of the model and fix boundary conditions and relations 
between parameters at unification energy. The constants /o, / 2 , Ja (momenta of 
the test function) are free parameters of the model: /q is related to the value of the 
coupling constants at unification, while the remaining two parameters / 2 and / 4 
enter in the expressions for the effective gravitational and cosmological constants 
of the model. The parameters o, b, c,5, e are functions of the Yukawa parameters 
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and Majorana mass terms of the model: 

a = Tr(YjY v + YjY e + 3(Y^Y U + Y^Y d )) 

b = Tr((YtY v f + (r e ty e ) 2 + 3(YjY u f + 3(yjY d ) 2 ) 

(2.4) c = Tr(MMt) 

= Tr((MMt) 2 ) 
c = Tt(MMtYjy„). 

3. RENORMALIZATION GROUP EQUATIONS IN NCG MODELS 

It is customary, within the various NCG models of particle physics coupled to 
gravity, to derive low energy estimates by assigning boundary conditions at unifi- 
cation compatible with the geometric constraints and running the renormalization 
group flow towards lower energies. The beta functions that give the RGE equations 
are imported from the relevant particle physics model. Thus, for exampls, the esti- 
mates derived in [5J (for the Connes-Lott model), or in [9] (for NC Yang-Mills), or 
in [28] and, more recently, in [13] are all based on the use of the RGE for the MSM. 
In the case of [13] an estimated correction term to the running for the top quark 
is introduced coming from the Yukawa coupling term for the r neutrino, which, in 
the presence of the Majorana mass terms, becomes of comparable magnitude. In 
[15] and [52], the same NCG model of [13] is analyzed using the full RGE analysis 
for the extension of the MSM by right handed neutrinos with Majorana masses, 
using the technique derived in [1], giving rise to different effective field theories in 
between the different see-saw scales. 



3.1. RGE equations for the i/MSM. The particle physics content of the CCM 
model :J3J is an extension of the minimal standard model (MSM) by right handed 
neutrinos with Majorana mass terms, which, as we recalled above, is a ^MSM 
model with unification energy. Renormalization group equations (at one loop) for 
this type of extension of the minimal standard model are described in pQ , while the 
renormalization group equations of the MSM are known at one and two loops [2]. 

The renormalization group equations for the r/MSM are a system of ordinary 
differential equations in unknown functions x — (xi) of the form 

(3.1) d tXi {t)=p Xi {x{i)), 

in the variable t — log(A/JWx), where the beta functions for the various parameters 
(coupling constants, Yukawa parameters, and Higgs quatric coupling) are given by 

41 o 19 , 7 , 

(3-2) * = 9ft?rf. ^-96^< A* = -16^*' 

for the three coupling constants, while the beta functions for the Yukawa parameters 
are of the form 

(3.3) 167T 2 p Yu = Y U ( 3 -Y^Y U - 3 -YjY d + a f Q g\ - \g\ - %g\) 

(3.4) 16^ 2 p Yd - Y d (^Y d - S -Y^Y U + a-\g\- \g\ - 8 ff 2 ) 

(3.5) 167T 2 p Yv = Y v {\y}Y v 3 2 Y}Y e + a - ~g\ - \&) 
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(3.6) 167T 2 fa = Y e (^Y e \YtY v + a~\gl- \&), 



with a as in (2.4 1 and where we use the notation 



~2 5 2 

The RGE for the Majorana mass terms has beta function 

(3.7) 16vr 2 p M = Y V Y}M + M{Y V Y}) T 
and the one for the Higgs self coupling A is given by 

(3.8) 16tt 2 & = 6A 2 - 3A(3. 9 | + g\) + 3gj +\{g\+ g\f + 4Ao - 8b 



where the terms a and b are as in (2.4 1 



Notice how, at one loop order, the equations for the three coupling constants 
decouple from the other parameters. This is no longer true at two loops in the 
MSM, see [2]. As explained in [1 , these equations can be treated as a series of 
equations for effective field theories between the different see-saw scales and inte- 
grated numerically. As shown in |18) , there is a sensitive dependence on the initial 
conditions at unification. 

3.2. Approximations of RGEs. We recall here some useful approximations to 
the RGEs above and some facts about the running of the solutions, with boundary 
conditions at unification compatible with the constraints of the CCM model recalled 
above. This analysis was already in [13] but we recall it here for direct comparison, 
later, with the modifications due to the interaction with the gravitational terms. 
The main idea is that, in first approximation, the top quark Yukawa parameter 



is the dominant term in the equations (3.3), (3.4), and (3.6). In the MSM case, it 



would also be dominant over all the terms in (3.5 1, but in the i/MSM, because these 



are coupled to the Majorana masses (3.7), the Yukawa coupling for the r neutrino 
can also be of comparable magnitude. 

If one neglects this additional term (that is, considers the RG equations of MSM) 
then in the running of the Higgs self-coupling parameter At one can also neglect all 
terms except the coupling constants gi and the Yukawa parameter of the top quark 
y t . This gives equations 

(3-9) & = ~ (24A 2 + 12Ay 2 - 9A( 5 2 + \gl) - 6 y 4 + ^ 2 4 + \g{ + \g\g\ 

where the Yukawa parameter for the top quark has 

(3-10) & = ^ - 8g£y - \g\y - y 2 g\y 

In [13] the problem of the existence of an additional term y T that contributes to 
the running is dealt with by modifying the boundary condition at unification for 
the top Yukawa parameter. This is dictated, in the NCG model, by the quadratic 
mass relation ( |6.3[ ), which gives 

£ {vZ? + (ve) 2 + z(y a u) 2 + m) 2 - 4. 9 2 . 

(regenerations 

Of these, the dominant terms become 

(3.11) z 2 + 3y 2 =4 5 2 , 
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Figure 1. Running of the coupling constants without anomalous 
dimensions (with "unification" at around 10 17 GeV). 



where Xt is the Yukawa coupling y° =3 (t) of the r neutrino (that is, for the generation 
index a — 3) and y t = yZ =3 (t) is the top Yukawa coupling. This gives a low 
evergy value yo for y t at yo ~ 1-04 as in [13] . computed assuming a value of the 
coupling constants at unification of g ~ 0.517 and setting the unification scale at 
A unif = 10 17 GeV. 



3.3. Coupling constants equations without gravitational corrections. The 

RGEs for the coupling constants (3.2) can be solved exactly. In fact, an ODE of 
the form 

(3.12) u'(t) = Au(tf, u(Q)=B, 
has exact nontrivial solutions 

(3.13) u(t) = ±— f= 

/ l-2AB 2 t 
V 

At A = Mz, the values of the coupling constants 

(3.14) 5l (0) = 0.3575, g 2 (0) = 0.6514, g 3 (0) = 1.221 

determine the values of B in the above equation, while A is determined by the beta 
functions of (3.2 1. This gives a running as illustrated in Figure [I] One obtains the 
usual picture where, at high energies, the coupling constants do not quite meet, 
but form a triangle. 



3.4. Top Yukawa coupling without gravitational corrections. With the so- 
lutions (3.13) to the running (3.2) of the coupling constants <7j(t), one can then 
solve numerically the equation for the running of the top Yukawa parameter, with 
beta function (3.10), using the boundary conditions at unification as in |13j . as 
described above. One obtains a running for y(t) as in Figure [2] The running of the 
Higgs quartic coupling can be analyzed in a similar way, by numerically solving the 
equation with beta function (3.10), see [T3] and [T5]. 
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Figure 2. Running of the top Yukawa coupling without anoma- 
lous dimensions. 



3.5. The running of the gravitational terms in the NCG models. The RGE 

analysis of the gravitational terms carried out in |13j is based on the running of the 
conformal gravity term derived in [3], while estimates are given on the parameter 
f 2 of the model, based on the value of the Newton constant at ordinary scales. In 
particular, the running of the conformal gravity term with the Weyl curvature 

«0 / Cuvpcr C^ f " T 



was analyzed in [T3] using as beta function (see equations (4.49) and (4.71) of [3]) 

a _ ^L133 2 

Pv ~ (4tt) 2 10 V ' 

with ao = It is shown in |13) that the running rj(t) does not change the value of 
T) by more than a single order of magnitude between unification scale and infrared 
energies. A different type of RGE analysis was done in [55], where the expressions 
(6.6) and (6.7| for effective gravitational and effective cosmological constants of the 



model were run, with the renormalization group equations of the i/MSM model, 
between unification and electroweak scales, to construct models of the very early 
universe with variable gravitational and cosmological constants. 

In the present paper, we propose a different use of the renormalization group 
analysis, based on introducing correction terms to the RGEs of the z/MSM model 
coming from the coupling of matter to gravity, through the method of Wetterich's 
functional renormalization group equations (FRGE) of [ID], as in [3"5] . 

4. Renormalization Group Equations and the Spectral Action 

As recalled in the previous section, generally in NCG models the RGE for the 
particle physics sector are imported from the corresponding particle physics model 
(MSM, as in [2], or an extension with right handed neutrinos and Majorana mass 
terms, see [T]). The fact that, in the NCG models, matter is coupled to gravity only 
enters through the constraints on the boundary conditions at unification described 
above, but not in the form of the RGEs themselves. We propose here a different 
RGE analysis, where an effect of the gravitational term is manifest also at the level 
of the RGEs. 
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4.1. Wetterich's FRGE and the spectral action. There is a general method, 
developed in [40] , to derive renormalization group equations from an effective action. 
In terms of the spectral action functional, one can regard asymptotic expansion for 
large A as the large cutoff scale limit of an effective action, somewhat like in [4"P] . 
Namely, we have a cutoff scale in the spectral action, if we think of the test function 
/ as a smooth approximation to a cutoff function. We can then interpret the 
spectral action Tr(/(D/A)) itself as the "flowing action" in the sense of Wetterich, 
with the dependence on A in f(D/A) moving the cutoff scale on D. In the limit 
for A — > oo, we can, to very high precision (see [TT]) replace the spectral action 
by its asymptotic expansion, which recovers the classical action for gravity coupled 
to matter. In the notation of the Wetterich approach we would then denote by 
Sk = Tr(/(_D/A)) the spectral action where f(x) is an even smooth approximation 
of a cutoff function with cutoff at x = k > 0. 

One could consider a functional renormalization group equation, in the sense of 
Wetterich, associated to the spectral action functional, in the form 

(-1.1! <V<V ■ = ;>( | r'^+Kk) W k ) , 




where t ~ log A and k is the cutoff scale for the test function /, the A = (Aj) are 
the boson fields, viewed as the independent fluctuations of the Dirac operator in 
the manifolds and non-commutative directions, and IZk is a tensorial cutoff as in 
[40) . A detailed discussion of how to derive RGEs from the spectral action with 
this method, the effect of adding the fermionic term of the spectral action as in 



(2.1 ) , and the relation between the RGEs obtained in this way to the RGEs of the 
z^MSM of [T| will be given elsewhere. 

For our purposes here, it suffices to observe that, for sufficiently large fc, the non- 
perturbative spectral action is very well approximated by an action of the from 

(4.2) 5[.g^]-^(^V(0 2 )-J-(0 2 )i?+^(^) 2 ) ^d 4 x, 

for suitable functions V and J- . Strictly speaking, this is true in the case of suf- 
ficiently regular (homogeneous and isotropic) geometries, for which the spectral 
action can be computed non-perturbatively, in terms of Poisson summation formu- 
lae applied to the spectrum of the Dirac operator, as shown in |llj . see also [7J, 
[23] . |24) . and [55] • In such cases, the functions V and T take the form of a slow-roll 
potential for the Higgs field, or for more general scalar fluctuations of the Dirac 
operator of the form D 2 H> D 2 + <fr 2 . These potentials have been studied as possible 
models of inflationary cosmology. A derivation of the non-perturbative spectral 
action and the slow-roll potential in terms of heat-kernel techniques is also given in 



For more general geometries, in adopting the form (4.2 1, as opposed to the more 



general (2.3), we are neglecting the additional gravitational term given by the Weyl 
curvature. While these can produce interesting effects that distinguish the NCG 
model of gravity from ordinary general relativity (see [35], [37J for a discussion of 
some interesting cases), the term is subdominant to the Einstein-Hilbert term at 
unification scale [52]. The argument from [T3J recalled above about the running of 
this term shows that it changes by at most an order of magnitude at lower scales, so 
we can assume that it remains subdominant and neglect it in first approximation. 
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4.2. RGE with gravitational corrections from the spectral action. The 

Wetterich method of FRGE was used in [T7] and to study the running of 
gravitational and cosmological constant for the Einstein-Hilbert action with a cos- 
mological term. The existence of an attractive UV fixed point was shown in [36] 
and further studied in [20]. This result was further generalized in [29] to the case 
of a theory with a scalar field non-minimally coupled to gravity via a pairing with 



the scalar curvature of the form (4.2) as above, where the values V(0) and ^(0) 



depend on the cosmological and gravitational constant, so that (4.2) also contains 
the usual Einstein-Hilbert and cosmological terms. The functions V and T are 
assumed to be polynomial in |29j . with V at most quartic and JF at most quadratic, 
although most of the results of [29] (see the generalization given in [25 ) work for 
more general real analytic V and J- ', which include the slow roll case. A detailed 
analysis of the critical surface of flow trajectories of the RGE approaching the UV 
fixed point is carried out in [29] . 



The explicit form of the RGE for an action of the form (4.2) were computed 
explicitly in [25], [29]. By the discussion above, they can serve as a good model of 
RGE for the spectral action functional. After making the same approximations as 



discussed in £3.2 the resulting equations become of the form used in [35], which we 



will discuss in detail in Sj5] below. 

4.3. RGE flow as a running of the noncommutative geometry. In terms 
of the geometry of the model, the renormalization group flow has a natural in- 
terpretation in terms of a running of the finite noncommutative geometry F = 
(Af,1-If, Dp). In fact, as shown in [13] (see also Chapter 1 of [IE]), the Dirac 
operator Dp of the finite geometry encodes all the Yukawa coupling matrices Y of 



(3.3)-(3.6), as well as the Majorana mass term matrix (3.7). Moreover, the data of 
the Dirac operator Dp on the finite geometry are described by a moduli space of 
the form M. p = Ge x Q q , with the quark sector given by 

G q = (C/(3) x CT(3))\(GL 3 (C) x GL 3 (C))/C/(3), 

and the lepton sector 

G e = (£7(3) x f/(3))\(GL 3 (C) x GL 3 (C) x 5 3 )/C/(3), 

with £> 3 the space of symmetric complex 3 x 3-matrices. (For a more general discus- 
sion of moduli spaces of finite spectral triples see also [BJ.) This space describes the 
bare parameters that enter the spectral action, in the form of the Dirac operator 
Da, with D — $x <8> 1 ffi 75 <8> Dp the Dirac operator of the product geometry and 
A the inner fluctuations. Thus, the one-loop RGE flow can be interpreted as a flow 
on A4p, for given solutions <7<(i) of the RGE for the coupling constants. 

This can be better interpreted as a flow on a fibration Ai p over a 3-dimensional 
base space Bp, with fiber the moduli space M.f, where the base space Bp cor- 
responds to the possible values of the coupling constants (31,52,93)- The renor- 
malization group flow then becomes a flow in M.p. At one-loop level, where the 
equation for the coupling constants uncouples from the equations for the Yukawa 
parameters and Majorana mass terms, one can see the solutions to the RGE flow 
for the coupling constants as defining a curve in the base space Bp and the RGE 
for the remaining parameters as defining a curve in M.p that covers the curve in 
Bp. 



ASYMPTOTIC SAFETY, HYPERGEOMETRIC FUNCTIONS AND HIGGS MASS 11 



5. RENORMALIZATION GROUP EQUATIONS WITH ANOMALOUS DIMENSIONS 

Asymptotic safety, as formulated in |39| . is a generalization of the notion of 
renormalizability, which consists of the requirement that the coupling constants lie 
on the critical surface of the UV fixed point. This allows for renormalization group 
analysis of gravitational terms coupled to matter. One of the main consequences 
of the asymptotic safety scenario is that the RGE for the matter sector acquires 
correction terms coming from the gravitational parameters in the general form 

(5.1) dtX] =^ M +(3f a \ 

?SM 



where Xj are the running parameters, t = log A, and /3j is the Standard Model 
Ej and /3| r 



beta function for Xj and /3| lav is the gravitational correction. The latter is of the 



for 



m 



a, A 2 



< 5 - 2 > v-^jm"- 

where the cij are the anomalous dimensions. 

The scale dependence of the Newton constant is given by (see [30] ) 

(5.3) Mj,(A) = M% + 2 Po A 2 , 

where the parameter p that expresses this scale dependence is estimated to have 
a value of p Q ~ 0.024, see [3DJ and [25J, [25]. 

We now consider the presence of the correction terms to the RGE coming from 



the gravitational sector, via the additional terms (5.2) in the beta functions (5.1) 



5.1. Anomalous dimensions in the gauge sector. As in [35], we assume that 
the cij for the couping constants are all equal and depending only on their common 
value at unification a± = 02 = 03 = a g , with \a g \ ~ 1 and a negative sign. Thus, 



instead of an equation of the form (3.12 ) for the running of the coupling constants, 



we obtain a modified equation of the form 
(5.4) u'{t) = -au(t) + Au(tf, 



u(0) = B, 



which again can be solved exactly and has (in addition to the trivial solution) 
solutions of the form 



(5.5) 



*(*) = 



I A + cxp 2« (I 



log(-^+^) 
2a 



Again imposing conditions (3.14) to fix B, and setting 

\aj 1 



167T/90 167T/90 



with A determined by the beta functions (3.2), with pa ~ 0.024 as above. We find 



that all the couplings run rapidly to zero at high energies, as illustrated in Figure 
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t log! 1(1) 



Figure 3. Running of the coupling constants with anomalous di- 
mensions. 



5.2. Anomalous dimension for the top Yukawa coupling. We have seen 
above that, at high energies, the coupling constant terms are then subdominant 
in the running for y and A for sufficiently large A. For the top Yukawa coupling, 
the analysis of 35j shows that the anomalous dimension satisfies a y < 0. Thus, one 
obtains an RGE equation for y t of the simpler form 



(5.6) 



otyt = -t^ — yt 



1 9 
16tt 2 2 



This can be solved exactly, with solutions that are again of the general form (5.5 ), 
this time with parameters 



(5.7) a = 

Again, the parameter is set to p 



167rpo 32-7H 

~ 0.024. There is an infrared fixed point with 
27r|a„| 



Vm 



9p 



We analyze the running of the top Yukawa coupling, in the presence of anomalous 
dimensions, using ( |5.6[ ) with the assumption, as discussed in [35] . that a g < a y < 0. 
The initial condition can be taken as in [T3|, with y(0) = B = 1.04. We obtain a 
running where y(t) rapidly decays to zero before reaching unification scale. One 
sees in this way the well known fact that the presence of the gravitational correction 
terms makes the Yukawa couplings asymptotically free. 



5.2.1. Anomalous dimension for the Higgs self- coupling. The qualitative analysis 
of the running of the Higgs self-coupling A given in |35j can be equally applied 
within the setting of the NCG model. They use an estimate on the anomalous 
dimension a\ ~ 3.1, as suggested in [55], [5S]. (We will discuss in § j]5.5| and 5.6 
some other natural choices of a\.) The effective action is of the form (4.2), which, 



as we discussed in f 4.2 above is a very good approximation for the spectral action 
for A running between the electroweak and the unification scale. The qualitative 
analysis of [5S] is based on estimating, as above, that the top Yukawa coupling 
contribution to the beta function for the Higgs self-coupling is dominant over the 
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tlog(lO) 



Figure 4. Running of the top Yukawa coupling with anomalous 
dimension \a v \ ~ 0.155. 



gauge contribution, which leads to a further simplification of the equation to 



(5.8) 



aX , . ) | \ 2 I • ) \ ,,2 ,, ' , 



16irp 



167T' 



-(24A 2 + I2\y 2 - 6y 4 



coupled to the running (5.6) of the top Yukawa y±- Here again po ~ 0.024. 



The equation (5.8) also gives an approximate fixed point for A determined by 
the equation 

Tra x 



(5.9) 
or equivalently, 
(5.10) 



24A 2 + l2Xy 2 IR 



-A = 0, 



Po 



7T ,8, . . , 47T 2 I a u 1 2 

24A 2 + -(-|a„| + a A )A - — = 0. 
Po •j 11 Po 



5.3. Riccati equations and singular linear ODEs. After substituting the ex- 



plicit general form of the solution for the equation (5.6) for the function y, the 



equation (5.8) for the unknown function A becomes a nonlinear ODE of Riccati 

type, 

(5.11) A' = qo{t)+ qi (t)\ + q 2 (t)\ 2 , 

in our case with 



(5.12) 



9o (t) = 



-3?/ 4 « 



<&(*) = 



3y 2 (t) 3 



8tt 2 ' ,iW 167rpo 4tt 2 ' 2tt 2 
Riccati equations have the property that they can be transformed into a linear 
second order ODE by a change of variables of the form 

v! 



(5.13) 

This gives a linear equation of the form 

«&(*) 
»(*) a = 



(5.14) ,/'_(g 1 (t) + 

Using a general form 



Ag 2 . 



u' + q 2 (t)q (t)u = 0. 



Ce 2at + A 
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for the solutions of (5.6 1, with parameters a, A, and C determined in terms of the 
coefficients of ( 5.6 1 and the initial condition, as discussed above, we obtain a second 

9 , a 



order linear equation of the form 
a\ 3 



(5.15) dtu{t) - ( 



167rp 4tt 2 Ce 2at + A 



)d t u 



16tt 4 v Ce 2at + A 



Yu(t) = 0. 



After changing variables to x = e 2at , with v(x) — v(e 2at ) = u(t), we can write 
the above as 
(5.16) 

a a 3 a w „ „ , NN 9 , a 



(2axydiv(x) + (2a- 



I6np„ Air 2 Cx + A 



)(2axd x v{x)) - 



1&tt^Cx + A' 



-v{x) = 0. 



5.4. Polynomial solutions of singular linear ODEs. We consider here special 
relations between the parameters a, C, A, and a\ that lead to the existence of 
special polynomial solutions of (5.16|. 

Following [31] , given a linear differential equation of the form 

(5.17) y" = X y' + s Q y, 

where Ao and so are C°° (I) functions on some interval of real numbers ICR, one 
can recursively obtain equations of a similar form for the higher derivatives, 

(5.18) y [n+1) = K-iy' + s n - lV , 

where the functions A„ and s n axe recursively defined by 

An = A__ 1 + s n -i + AoA„_i 

(5.19) _ ; 



Then the equation (5.171 has a polynomial solution of degree at most n if A n A„_i ^ 
and 



(5.20) S n := A„ s n _i - A„_i s T . 

One then has a solution of the form 



(5.21) 



y = exp 



exp 



0. 



gn-l 

A n -i 



and a general solution of the form 



(5.22) 



!J = exp 



C x + C 2 / exp / A + 2 



'A,, 



The equation (5.161 can be put in the form (|5.17|) with 
(5.23) Ao 



3a 
"4tt 2 



2a- 



16-Kpo 



(A + Cx) 



(5.24) 



so = 



2ax(A + Cx) 
9 

'6A7T 2 x 2 (A + Cxy 
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The first step of the recursive equation (5.191 gives 
(5.25) 

x - i (r 2 \ ° 2a2 >- .. c2a * \ 

Al ~ (A+Cx) 2 ~r 1024a 2 7r 2 p 2 Wa-n Po ) 



1 f a2 , 9 3A | A 2 a| , 3Aa A _ A 2 a A 9_' 

(A+Ci) 2 i 2 \^ T 64-Ti- 4 4tt 2 ' 1024a 2 7r 2 p 2 ' 128a7r 3 p 16a7rp 64tt 2 

1 f A I 3 I Aoa 

(A+Cx)x 2 \ ^ ~ 8tt 2 t 32airp 



(5.26) 



(A+Cx) 2 x 4tt 2 ^ 512a 2 7r 2 p 2 ^ 128a7r 3 p 

. . 27 | 9A , 9Aa A 

1 (A+Ca;) 3 a; 3 \^512tt 4 64tt 2 2048a7r 3 p 

1 / 27C , 9Ca x \ 

T (A+Cx^x 2 \64-x 2 2048a7r 3 p 0/ l • 

Thus, at this first step one obtains 

1 / 81A 27A 2 9A 3 9A 3 a x \ 

dl ~ x A {A + Cxf V40967T 4 + 512tt 4 ~~ 64tt 2 + 2048arr 3 po J 

1 / 81C 21 AC 27 A 2 C 27A 2 Ca x \ 

x 3 (A + Cx) 5 V40967T 4 + 128tt 4 ~ 32tt 2 + 1024a7r^ / + 

1 / 81C 2 81AC 2 27 AC 2 a x \ 

x 2 {A + Cx) 5 V5127T 4 _ 64rr 2 + 1024a7r 3 p / + 

1 / 9C* 3 9C 3 a A 



x(A + C*a;) 5 V 16vr 2 512a7r 3 p , 
The vanishing <5i = implies that the parameters should satisfy either C = and 
A = or else C = and 

(-9 - 24A + 64A 2 ^ 2 )qp 
(5.27) a A = — . 

Since both these solutions have C = they are ruled out for physical reasons. We 
nonetheless look more closely at the second case. The linear solution in this case 
is obtained by replacing these values of C and a\ in (5.23) and (5.24), so that one 
gets 

so | 1 

\ r-n n _ (-g- 24 ^ + 64A 2 ^ 2 )qp — 

A solution is then given by y = exp(J Sq/Xq) = x and the general solution by 

y = Cix + C 2 x a+2 , 

where 

9 

a ~ + 64A^' 
with the latter obtained from the fact that 

9 

Ao 



„_ _ (-9-21_,H-01.V7H)ap — 2 • 

G-0,a A - 64:A z Tr z X 

At the second order, the explicit expressions for A2 and s 2 are more involved, 
but one encounters essentially the same situation. Namely, the condition for the 
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vanishing 82 = is given by a system of five equations in the parameters. These 
equations have solutions 

(-9- 24A + 64A 2 7r 2 )ap 



C = 0, 



ax 



C = 0, 



2A 2 TT 

3ap (3+ 16A) 

4A 2 7T ' 



which again are not compatible with having C ^ in the solution of the RGE 
equation for the top Yukawa coupling y. 

One can see a similar situation for the case of A3 and S3, where one has several 
more polynomial solutions, but all of them again under the condition that C = 0. 
In this case one obtains a system of seven equations for the parameters that arise 
from imposing £3 = 0. The possible solutions to these equations with A positive 
real are the following: 

a (-9 - 24A + 64A 2 n 2 ) r 

ax= 2A^r 

a (3 + 2AA + 64^ 2 tt 2 ) r 
a\ = — 



(5.28) 



.4 



A = 



24 + V576 + 4608tt 2 



48- 



256tt 2 
^2304 + 



2A 2 TT 

3ap (3 + 16A) 

4A 2 7T 

-32a7rp 



C = 



4608tt 2 



a\ = — 32a7rpo 



C = 



C = 



C = 



256tt 2 

Notice how, at least up to degree three, the only condition that is compatible 



with having a positive anomalous dimension a\ > is the relation (5.27) (though 
not for all values of A). 

5.5. Higgs self-coupling and hypergeometric functions. We have seen above 
that, in the non-physical case with C — and with the anomalous dimension a\ of 



the form (5.27), the equation (5.16) admits linear solutions y = C\x. Here we drop 



the unphysical condition C = 0, but we keep the same constraint on a\ given (5.27), 



and we show that then the resulting equation (5.16) can be integrated in closed 



form and has a general solution that can be expressed in terms of hypergeometric 
functions. 



Notice first that, after imposing the condition (5.27), the equation (5.16) can be 



written in the form (5.17), with 

1 



(5.29) 



An 



so 



9C 



(A 



Cx) V64A 2 tt 2 
9 



3C 
8Att 2 



64An 2 x 



64ir 2 x 2 (A + Cx) 2 
We introduce the following auxiliary functions: 
(5.30) h{x) = (A + Cx) a ~ 3 

with 

, ^9 



(5.31) 



4vr 2 (-9 + 4A(3 + 4A7r 2 )) 
WAn 2 
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-2 



Figure 5. Solutions of the form (5.371 with v(x) as in (5.36) can 
exhibit complicated and singular behavior, depending on the values 
of the parameters C\, C 2 , A, C, a. 



and 
(5.32) 

(5.33) 

We also set 
(5.34) 

(5.35) 



16Att 2 



Cx 



with 77 = 



64A 2 7r 2 ' 



p = a — r) — 1 



8An 2 



Cx 
~A' 



Then the equation (5.17) with (5.29) has general solution of the form 
v{x)= Cxh^htx) 2 F 1 ( a ,/3,2- V ,~,(x)) 

(5.36) 

+ C 2 x*h{x) 2 F 1 (a-(2+^),p + (2+^)r } ,r 1 ^(x)), 
where 2 Fi(a, b, c, z) is the Gauss hypergeometric function. 



The corresponding solutions for the original Riccati equation (5.8) that gives the 



renormalization group flow of the Higgs self coupling A are then of the form 

2tt 2 «'(t) 



(5.37) 



A(t) 



with u{t) = v(e 2at ) 



3 u(t) 

As was already observed in |35| , these solutions, for varying choices of the param- 
eters Ci, C2, A, C, a, can exhibit very complicated and singular behavior (see an 
example in Figure [5]). 

Thus, one needs to ensure the existence of particular solutions that will be com- 
patible with physical assumptions on A. So at least one needs a solution for which 
A(i) is positive and does not develop singularities for all t > 0. This is, in fact, easy 
to achieve. For example, setting A = 1 and C > 2, and a = \a y \ / (I6ir po) , with 
\a y \ = 0.155 and po = 0.024, gives a solution of the form as in Figure [6j 

However, solutions obtained in this way are still not entirely satisfactory from 
the physical point of view for two main reasons. The first is that we are using the 



18 



CHRISTOPHER ESTRADA AND MATILDE MARCOLLI 




Figure 6. A solution with X(t) continuous and positive for t > 0. 



relation (5.27) that expresses a\ as a function of A, a and po- While the example 
A = 1 above gives a plausible value of the anomalous dimension a\, the value 



A = considered in (5.7) would give a\ the wrong sign. The other reason is 



that a solution like the one illustrated in Figure [6j while well behaved in the range 
t > 0, is still plagued by infrared divergences, although these can be pushed very 
far back in the t < range. 

5.6. Other special solutions of the Higgs self-coupling ODE. The form of 



the singular linear ODE (5.16) also suggests as a possible natural choice of a\ the 
value 



a\ = 'd2np a a = 2 \a y 



(5.38) 

For this value the equation (5.16) simplifies to the form 

3 a 



(5.39) {2axyd 2 x v{x) + 



4ir 2 Cx + A 



2axd x v{x) 



9 



16tt 4 v Cx + A 



fv{x) = 0, 



or equivalently to an equation of the form (5.17), with 

3 



(5.40) 



An= - 



•so 



8?r 2 a:(A + Cx) 
9 



2 \ 2 



64ir 2 x 2 (A + Cx) 2 ~ " 

In this case, again, we have an explicit family of solutions in terms of hypergeometric 
functions. We introduce the following notation. 

-20 



(5.41) 

with 
(5.42) 



h x {x) = {x{A + Cx)f, h 2 (x) = 



Cx 
~A 



l6An 2 



and 9 as in (5.31 1 and 
(5.43) 6 



y/9 + 4tt 2 (-9 + 4A (-3 + AAtt 2 )) 
16Ait 2 
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Figure 7. A solution v(x) with a\ = 32irpoa and with A = 
9/(327r 2 ) and C — 2, with oscillatory behavior. 



(5.44) 



We then have a general solution of the form 
v(x)= Cih x {x) 2 F 1 {0 + 9,1 + 
+ C 2 h 1 (x)h 2 {x) 2 F 1 {6- 



+ M + 20,7(*)) 

i + e-e,i-2e,j(x)), 



with 7(2;) as in (5.351, which is again given in terms of Gauss hypergeometric 
functions. 

One advantage with respect to the family of solutions discussed in the previous 
section is that now the value A 



32ir 



1 of (|5.7|) is compatible with a positive 



anomalous dimension a\ . In this case again one can construct examples of solutions 
u(t) = v(e 2at ) for which the resulting function \{t) as in (5.37) is continuous and 
positive for all t > 0. 



5.7. Implementing a Shaposhnikov Wetterich-type Higgs mass estimate. 

The Higgs mass estimate of 126 GeV derived by Shaposhnikov Wetterich in [55] . 



uses the RGE with anomalous dimensions (5.6) and (5.8), with values of the pa- 
rameters A = 9/(327r 2 ), a = |aj,|/(167rpo) with p ~ 0.024 and \a y \ ~ 0.155, 
C = -gz — A, with B = y(0), where the value y(0) is constrained physically by 
compatibility with the top quark mass. In [35] they use a value of a\ ~ 3.1 for 
the anomalous dimension, though they note that their qualitative argument does 
not depend on the precise values of the anomalous dimensions. Here we work with 
the explicit solutions in terms of hypergeometric functions, which are instead very 
sensitive to the explicit values of a\ and a y . We show here that a slightly larger 
value of a\ ~ 5.08 is compatible with the value of \a y \ ~ 0.155 used in [35] and 
with the constraint on y(0) imposed by compatibility with the top quark mass. 
In particular, for rrit ~ 171.3 GeV, and the Higgs vacuum vjj ~ 246 GeV, which 
means a value y(0) ~ 0.98, and keeping the same value of A = 9/(327r 2 ), we obtain 
C - 0.09. 



We have a solution of the equation (5.8) of the form (5.37), with u(t) = v(e 2at ) 



and v(x) a solution of (5.16) of the form 
/ 9 



(5.45) v(x) = %t(x° 



V32C7T 2 



32C7rV 



2/3 



-feV*)), 
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Figure 8. A solution A (a;) with a A = 5.08 and with A = 9/(327r 2 ), 
compatible with a top quark mass of m t ~ 171.3 GeV and giving 
a Higgs mass to# ~ 125.4. 



for 



a = - 39) + \ 1 3 + J- 126 + 96 2 + (144 - 2304tt 2 ) 
6 \ V 36 



uj= ~(1 + ~(-3 + 36 + \/4 - 126 + 96 2 - 64tt 2 - 2-^-49 + 64vr 2 ) 



i(-3 - W-126 + 96 2 + —(144 - 2304tt 2 )) 
3 V 36 



12 



i(18 - 66 - V-196 + 2567T 2 + 2W-126 + 96 2 + —(144 - 2304tt 2 ))) 



36 



- f-3 + 36 + \A - 126 + 96 2 - 64tt 2 - i ^-49 + QAn 2 ) 



if 18 - 66 - iV-196 + 256tt 2 + 2 J -126 + 96 2 + — (144 - 2304tt 2 ) | 
12 \ V 36 J 




126 + 96 2 + — (144 - 2304?r 2 ) . 

36 I 

Here the parameter A has been assigned the value A = 9/(327r 2 ), while, as above, 
the parameter C is given by C = — A, with B — y(0) and the parameter 6 
is given by imposing a relation of the form a\ — 327rpoa6. With the parameters 
assigned as above, we find that A(i) is continuous and positive for alH > and such 
that A(0) = 0.130, compatibly with a Higgs mass estimate of m# = V2Xvh ~ 125.4, 
for a slightly larger anomalous dimension a\ ~ 5.08. The behavior of this solution 
for x > 1 (that is, for t > 0) is illustrated in Figure |8j where we plotted the function 

47r 2 a v'(x) 



X(x) 



v(x) 



ASYMPTOTIC SAFETY, HYPERGEOMETRIC FUNCTIONS AND HIGGS MASS 21 



5.8. Infrared behavior. The solutions obtained in this way acquire non-physical 
negative values, and possibly divergences, at the infrared end of the spectrum. This 
phenomenon was noticed also in |35j and related to radiatively induced spontaneous 
symmetry breaking of Coleman- Weinberg type [14] . In fact, the problem with the 
behavior of solutions for large negative t < can simply be attributed to the 
failure, in this range, of the approximations (5.6) and (5.8), where the terms with 



the coupling constants should be reintroduced as in (3.9) and (3.10), with explicit 
solutions for the coupling constants of the form (5.5). 



6. Geometric conditions at unification 

In the models of particle physics based on noncommutative geometry and the 
spectral action functional, there are several constraints on the parameters at unifi- 
cation energy, imposed by the geometry of the model. We recall here the constraints 
derived within the model of |13) and we discuss how these are affected by the pres- 
ence of the gravitational correction terms discussed above in the renormalization 
group flow. 



6.1. Constraints at unification in the CCM model. As discussed in [13] and 
in §1 of |16) . there are constraints on the values of the parameters listed above at 
unification energy, imposed by the geometry of the model. We recall them here and 
we refer the reader to the references above for a detailed account of how one obtains 
these relations. We then show that, in the presence of anomalous dimensions the 
geometric conditions at unification can be expressed in terms of the presence of a 
Gaussian matter fixed point Fjjv — (Af,"W.f,0) for the flow on the moduli space 
M.F of the finite geometries, and this changes the interpretation of the constraints 
at unification energy with respect to [13] and [16]. In particular, the Majorana 
mass terms become compatible with the dark matter models of Shaposhnikov and 
Tkachev [33] , [33J , [33] , but one loses the geometric see-saw mechanism of [T3J . 



6.1.1. Coupling constants. The normalization of the Yang-Mills terms in the as- 
ymptotic expansion of the spectral action ([13]) leads to a unification of the coupling 
constants with the common value g at unification energy satisfying 

This fixes the value of the parameter fo of the model in terms of the value at 
unification of the coupling constants g 2 = gf|A unif = fflUunit = (5/3)Si |A unif ■ 



6.1.2. Higgs vacuum. The parameter a is related to the Higgs vacuum vh, through 
the relation 



(6.2) =v H , 

which follows from a direct comparison of the terms involving the Yukawa coupling 
in the asymptotic expansion of the spectral action and in the Lagrangian of the 
Standard Model (see (4.13) and (4.15) of [13]). 
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6.1.3. Mass relation. The two relations above are used in 13J to derive an addi- 
tional quadratic mass relation at unification, which follows again from (4.13) and 
(4.15) and the expression vjj — 2Mw/g for the Higgs vacuum, where M\y the 



W-boson mass. One obtains then, using (6.1| and (6.2), 

(6.3) ^ ( m i + m e + 3m « + 3TO d)|A=A u „if = 8Af wU=A uni f 

generations 

6.1.4. Higgs self-coupling. The value of the Higgs self-coupling parameter A at uni- 
fication energy is predicted to satisfy the relation 

7T 2 b(A U nif) 



(6.4) A(A U 



2/oa(A unif ) 2 - 



6.1.5. See-saw mechanism. The see-saw mechanism for the neutrino sector, imple- 
mented geometrically by the component of the Dirac operator on F acting on the 
right-handed neutrino sector of the fermion representation, imposes a condition at 
unification on the parameters c, fa and /o 

(6.5) < c(A unif ) < 6 J^. 

To Jo 

6.1.6. Gravitational constant. At unification energy the (effective) gravitational 
constant is related to the parameters c, f% and /o of the model through 

37T 

(6 ' 6) Gcff(Aunif) = 192/ 2 A? mif -2/oc(A unif )- 

6.1.7. Cosmological constant. At unification energy the (effective) gravitational 
constant is related to the parameters of the model through 

(6.7) 7eff(A unif ) = (192 h Knit ~ 4/ 2 AL f c(A unif ) + /„ D(A unif )). 

(We use here the notation 7 for the cosmological constant, as the more commonly 
used A is already taken by the standard notation in NCG for the energy scale in 
the spectral action.) 

6.1.8. Conformal gravity terms. The conformal gravity terms given by the Weyl 
curvature term J C illlpa C llv( " J and the conformal coupling of the Higgs field to 
gravity, appear, respectively, with coefficients a = — 3/o/(107r 2 ) and £0 = 1/12. 
The latter is the standard coefficient for the conformal coupling of the Higgs to 
gravity. These coefficients are computed after a normalization of the kinetic term 
of the Higgs. 

6.2. The case with anomalous dimensions. When we consider renormalization 
group equations with gravitational corrections given by anomalous dimensions as 
in f|5j these have a Gaussian matter fixed point at high energies, where the matter 
fields becomes uncoupled and asymptotically free. This means that the coupling 
constants and the Yukawa parameters tend rapidly to zero at high energies, as 
we have seen in the behavior of the solutions to the approximate RGE equations 
discussed in <|5j This leads to a corresponding reinterpretation of the geometric 
boundary conditions at unification. 

In fact, if we interpret the RGE flow as a flow of the finite geometry F, the 
Gaussian matter fixed point corresponds to a fixed point in the moduli space Mp 
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as described in §4.3| This fixed point has trivial coupling constants, and in the 
fiber Aip over this point of the base Bp it also has trivial Yukawa couplings, so 
that Fjjv corresponds to the finite geometry Fjjy = (Af, T~Lf, 0) with trivial Dirac 
operator. Notice that D — is also a fixed point of all the inner fluctuations of the 
geometry. Thus, the geometric constraints at unification now become the condition 
that the finite geometry is close to the fixed point Fjjv ■ 

In terms of physical consequences of these boundary conditions, one of the pri- 
mary effects is on the Majorana mass terms. In the setting of [13] the top eigenvalue 
of the Majorana mass terms matrix M is assumed to be near unification scale, so 
that the three see-saw scales fall in between the unification and the electroweak scale 
(see a detailed discussion of this point in This is convenient for a geometric 

interpretation of the see-saw mechanism for neutrino masses (see [13] ). In the set- 
ting with anomalous dimensions and Gaussian matter fixed point the behavior of 
the Majorana mass terms becomes instead consistent with the dark matter models 
of Shaposhnikov and Tkachev [32], [33], [31], where the Majorana masses are at or 
below the electroweak scale, thus making them plausible dark matter candidates, 
in agreement with the exclusion estimates of Kusenko [19 . 

Comparing with the boundary conditions of |13j recalled above, we see that in 



the case with Gaussian matter fixed point F\jy = („4f,%f,0), the relation (6.4) 



can be used to fix /o, while (6.2) fixes the shape of the quartic potential near F\jy. 
The relation (6.5) still fixes /2, but now no longer implies a see-saw mechanism 
for neutrino masses, for the reason mentioned above that the Majorana mass scales 
can fall below the electroweak scale as in the Shaposhnikov-Tkachev models of dark 
matter [35], [33J, [31]. 

7. Further questions 

Some interesting questions arise in connection to the RGE equations discussed in 
the previous sections. It would be desirable to have a full derivation of Wetterich's 



FRGE (4.1 ), directly expressed in terms of the non-perturbative form of the spectral 



action, and their compatibility with the RGE flow of [T] for the particle physics 



sector of the model, with gravitational corrections as in (5.1), (5.2). In view of 
the recent interest in cosmological applications of the NCG models, it would also 
be interesting to see how the gravitational corrections to the RGE flow affect the 
RGE analysis of early universe models of [55] and the constraints on the use of 
the Higgs as a possible inflation mechanism within this type of models, see [5]. A 
related question is whether one can obtain constraints on the cosmic topology, via 
the analysis of [53] , [51] , and [7] , and the effects of the gravitational terms on the 
running. Other NCG models have been proposed recently, among them the model 
with supersymmetry of [1] and a model with one lighter and one heavier Higgs |37j . 
It would be interesting to see what effect the inclusion of gravitational corrections 



of the form (5.1 1, (5.2) has on the RGE analysis for these models. In particular, the 
supersymmetric model of @], first developed only for the QCD sector, will likely 
extend to the electroweak sector in a way that is not identical to the MSSM, hence 
the RGE equations will likely be also different from those of MSSM. The Wetterich 
method may provide a useful way to perform the RGE analysis of the spectral 
action for that model. 



24 



CHRISTOPHER ESTRADA AND MATILDE MARCOLLI 



Acknowledgments. The first author was supported by a Summer Undergraduate 
Research Fellowship. The second author is partially supported by NSF grants DMS- 
0901221, DMS-1007207, DMS-1201512, PHY-1205440. 



References 

[1] S. Antusch, J. Kersten, M. Lindner, M. Ratz, M.A. Schmidt Running neutrino mass param- 
eters in see-saw scenarios, JHEP 03 (2005) 024, |hep-ph/0501272 y3. 

[2] H. Arason, D.J. Castano, B. Kesthlyi, E.J. Piard, P. Ramond, B.D. Wright, Renormalization- 
group study of the standard model and its extensions: the standard model, Phys. Rev. D, 46 
(1992) N.9, 3945-3965. 

[3] I.G. Avramidi, Covariant methods for the calculation of the effective action in quantum 

field theory and investigation of higher- derivative quantum gravity, PhD Thesis, Moscow 

University, 1986, |hep-th/9510140| 
[4] T. van den Broek, W.D. van Suijlekom, Supersymmetric QCD and noncommutative geometry, 

Comm. Math. Phys. 303 (2011), no. 1, 149-173. 
[5] M. Buck, M. Fairbairn, M. Sakellariadou, Inflation in models with conformally coupled scalar 

fields: An application to the noncommutative spectral action, Phys. Rev. D 82 (2010) 043509 

[14 pages] 

[6] B. Cacic, Moduli spaces of Dirac operators for finite spectral triples, in "Quantum groups 
and noncommutative spaces", pp. 9—68, Aspects Math., E41, Vieweg Teubner, Wiesbaden, 
2011. 

[7] B. Cacic, M. Marcolli, K. Teh, Coupling of gravity to matter, spectral action and cosmic 

topology, arXiv:1106.5473 to appear in Journal of Noncommutative Geometry. 
[8] L. Carminati, B. Iochum, T. Schucker, The noncommutative constraints on the standard 

model a la Cannes, J. Math. Phys. 38 (1997), no. 3, 1269-1280. 
[9] L. Carminati, B. Iochum, D. Kastler, T. Schucker, Relativity, noncommutative geometry, 

renormalization and particle physics. Coherent states, differential and quantum geometry, 

Rep. Math. Phys. 43 (1999), no. 1-2, 53-71. 
[10] A. Chamseddine, A. Connes, The spectral action principle, Commun. Math. Phys. 186 (1997) 

731-750. 

[11] A. Chamseddine, A. Connes, The uncanny precision of the spectral action, Comm. Math. 
Phys. 293 (2010), no. 3, 867-897. 

[12] A. Chamseddine, A. Connes, Resilience of the Spectral Standard Model, arXiv:1208.1030 

[13] A. Chamseddine, A. Connes, M. Marcolli, Gravity and the standard model with neutrino 
mixing, Adv. Theor. Math. Phys. 11 (2007), no. 6, 991-1089. 

[14] S. Coleman, E. Weinberg, Radiative corrections as the origin of spontaneous symmetry break- 
ing, Phys. Rev. D7 (1973) 1888-1910. 

[15] A. Connes, Gravity coupled with matter and the foundation of noncommutative geometry, 
Comm. Math. Phys. 182 (1996), no. 1, 155-176. 

[16] A. Connes, M. Marcolli, Noncommutative geometry, quantum fields and motives, Colloquium 
Publication, Vol.55, American Math. Soc. 2008. 

[17] D. Dou, R. Percacci, The Running Gravitational Couplings, Class. Quant. Grav. 15 (1998) 
3449-3468. 

[18] D. Kolodrubetz, M. Marcolli, Boundary conditions of the RGE flow in the noncommutative 
geometry approach to particle physics and cosmology, Physics Letters B 693 (2010) 166-174. 

[19] A. Kusenko, Sterile neutrinos: the dark side of the light fermions, arXiv:0906.2968 

[20] O. Lauscher, M. Reuter, Ultraviolet fixed point and generalized flow equation of quantum 
gravity, Phys. Rev. D 65 (2002) 025013 [44 pages] 

[21] M. Marcolli, Building cosmological models via noncommutative geometry, International Jour- 
nal of Geometric Methods in Modern Physics, Vol.8, N.5 (2011) 1131-1168. 

[22] M. Marcolli, E. Pierpaoli, Early universe models from noncommutative geometry, Adv. 
Theor. Math. Phys. 14 (2010), no. 5, 1373-1432. 

[23] M. Marcolli, E. Pierpaoli, K. Teh, The spectral action and cosmic topology, Comm. Math. 
Phys. 304 (2011), no. 1, 125-174. 

[24] M. Marcolli, E. Pierpaoli, K. Teh, The coupling of topology and inflation in noncommutative 
cosmology, Comm. Math. Phys. 309 (2012), no. 2, 341-369. 



ASYMPTOTIC SAFETY, HYPERGEOMETRIC FUNCTIONS AND HIGGS MASS 25 



[25] G. Narain, R. Percacci, Renormalization group flow in scalar-tensor theories, I, 
larXiv:0911.0386t r2. 

[26] W. Nelson, J. Ochoa, M. Sakellariadou, Constraining the noncommutative spectral action via 
astrophysical observations, Phys. Rev. Lett. 105 (2010) 101602 [5 pages] 

[27] W. Nelson, J. Ochoa, M. Sakellariadou, Gravitational waves in the spectral action of non- 
commutative geometry, Phys. Rev. D 82 (2010) 085021 [15 pages] 

[28] Y. Okumura, Renormalization group analysis of the Higgs boson mass in a noncommutative 
differential geometry, Prog. Theor. Phys. 98 (1997) 1333-1344. 

[29] R. Percacci, D. Perini, Asymptotic safety of gravity coupled to matter, hep-th/0304222 

[30] M. Reuter, Nonperturbative evolution equation for quantum gravity, Phys Rev D 57 (1998) 
N.2, 971-985. 

[31] N. Saad, R.L. Hall, H. Ciftci, Criterion for polynomial solutions to a class of linear differential 

equations of second order, J. Phys. A: Math. Gen. 39 (2006) 13445-13454. 
[32] M. Shaposhnikov, The uMSM, dark matter and baryon asymmetry of the Universe, Journal 

of Physics: Conference Series 39 (2006) 9-11. 
[33] M. Shaposhnikov, The vMSM, leptonic asymmetries, and properties of singlet fermions, 

JHEP (2008) 0808:008, 56pp. 
[34] M. Shaposhnikov, I. Tkachev, The vMSM, inflation, and dark matter, Phys. Lett. B Vol.639 

(2006) 414-417. 

[35] M. Shaposhnikov, C. Wetterich, Asymptotic safety of gravity and the Higgs boson mass, 
larXiv:0912.0208t r2. 

[36] W. Souma, Non-trivial ultraviolet fixed point in quantum gravity, Prog. Theor. Phys. 102 
(1999) 181-195. 

[37] C.A. Stephan, New scalar fields in noncommutative geometry, Phys. Rev. D 79 (2009) 065013 
[38] K. Teh, Nonperturbative Spectral Action of Round Coset Spaces of 577(2), arXiv:1010.1827| 

to appear in Journal of Noncommutative Geometry. 
[39] S. Weinberg, Ultraviolet divergences in quantum theories of gravitation, in "General Relativ- 
ity: an Einstein centenary survey" (S.W. Hawking, W. Israel, Eds.) Cambridge Univ. Press, 
1979, pp. 790-831. 

[40] C. Wetterich, Exact evolution equation for the effective potential, Phys. Lett. B 301 (1993) 
90-94. 

[41] O. Zanusso, L. Zambelli, G.P. Vacca, R. Percacci, Gravitational corrections to Yukawa sys- 
tems, |arXiv:090 4.0938v2. 

Division of Physics, Mathematics and Astronomy, Mail Code 253-37, Caltech, 1200 
E. California Blvd. Pasadena, CA 91125, USA 
E-mail address: c.estrada@caltech.edu 
E-mail address: matildeScaltech.edu 



